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1 Introduction 

It is well known the meaning of methods of the functional differentiation and functional inte- 
gration for the problems of mathematical physics [HE] and applications [3j, in particular with 
regard to many-particle systems of statistical mechanics or quantum field theory. As far as we 
know for the first time equations in functional derivatives for the description of the evolution of 
states of many-particle systems were used by T.H. Gronwall [2], M.A. Leontovich [5] and N.N. 
Bogolyubov [6j[7]. Bogolyubov's results were developed by R.L. Lewis in the paper pj, where it 
was constructed the reduced representation of a nonperturbative solution of the Cauchy prob- 
lem of the BBGKY hierarchy. On quantum many-particle systems Bogolyubov approach was 
extended in [9]. 

In this article we present the rigorous results on the evolution equations in functional deriva- 
tives for generating functionals of states and observables of classical many-particle systems, 
namely the BBGKY hierarchy and the dual BBGKY hierarchy in functional derivatives respec- 
tively [in]. One more approach to the description of the evolution of states of many-particle 
systems is given by means of generating functionals for correlation functions. The generating 
functionals of correlation functions and marginal correlation functions are governed by the von 
Neumann hierarchy and the nonlinear BBGKY hierarchy in functional derivatives respectively. 
On the basis of the developed approach we construct nonperturbative solutions of the Cauchy 
problem of corresponding evolution equations in appropriate Banach spaces. 
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We briefly outline the results and structure of the paper. In section 2 we set forth an 
approach to the description of the evolution of states of many-particle systems within the 
framework of the generating functional of marginal distribution functions governed by the 
BBGKY hierarchy in functional derivatives. Then by means of introduced generating functional 
of marginal distribution functions we construct a nonperturbative solution of the initial-value 
problem of the BBGKY hierarchy on the basis of cluster expansions of groups of operators 
which describe dynamics of finitely many particles. 

In section 3 we develop an alternative method to the description of the evolution of states 
within the framework of the generating functionals of correlation functions and marginal corre- 
lation functions governed by the von Neumann hierarchy and the nonlinear BBGKY hierarchy 
in functional derivatives respectively. The relationships of the introduced generating functionals 
of states are established. 

In section 4 we state an equivalent approach to the description of the evolution of many- 
particle systems within the framework of the generating functional of marginal observables 
governed by the dual BBGKY hierarchy in functional derivatives. A nonperturbative solution 
of the initial- value problem of the dual BBGKY hierarchy is constructed on the basis of cluster 
expansions of groups of operators which describe the evolution of observables. 

Finally in section 5 we conclude with some observations and perspectives for future research. 

2 The BBGKY hierarchy in functional derivatives 

2.1 A generating functional of marginal distribution functions 

Let / = (/o, fi(xx), . . . , f n (%i, ■ ■ ■ , x n ), . . .) be a sequence of measurable functions defined on 
the phase spaces of the corresponding number of particles and u = (1, u(xi), . . . , YYi=i u i x i)i ■ ■ •) 
is a sequence of products of continuously different iable integrable functions u(xi), i > 1. The 
generating functional of functions / n , % > 1 is defined by the following expression 

(/,«) = 5Z - ! / •" / fn(xi,...,x n ) Y[u(xi)dxi . . . dx n . (1) 

For the integrable functions u(xi), i > 1, functional (JT]) exists. We refer to functional ([T]) as 
the generating functional of functions f n (xx, . . . , x n ), n > 1, by reason of the equality 

f n (x 1 ,...,X n ) = — r{f,U) \ u=0 , (2) 

du(xi) . . . ou[x n ) 

where S n /Su(xi) . . . Su(x n ) is the nth order functional derivative (the nth order Gateaux deriva- 
tive puifa]). 

We consider a system of a non-fixed number of particles in the space R 3 defined by the 
Hamiltonian [T3] 

n n 

H n = Y,K(Vi)+ E *(*-«&)> ( 3 ) 
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where K(pi) is a kinetic energy of the % particle and $(g«— qj) is a two-body interaction potential. 
Hereafter we assume that the function $ satisfies conditions which guarantee the existence of a 
global in time solution of the Hamilton equations for finitely many particles [13]. The evolution 
of states of a system is described by the sequence D{t) = (1, Di(t, Xi), . . . , D n (t, x±, . . . , x n ), . . .) 
of distribution functions. The integrable function D n (t,xi, ■ ■ ■ ,x n ) is defined on the n-particle 
phase space and symmetric with respect to the permutations of the arguments xi, . . . , x s . We 
introduce the generating functional (D(t),u) of the distribution functions D n (t, x±, . . . , x n ), n > 
1, according to definition (pQ) by the expansion 

oo _ „ „ n 

(D(t),u) = ^2~ ■■■ D n (t,xx,...,x n )Y[u(xi)dxi...dx n . (4) 

n=0 U - J J i=l 

The evolution of states of many-particle systems in terms of generating functional (jlj) is 
governed by the following evolution equation in functional derivatives 

j t (D(t),u) = J {Kip,), *_(£(«), «)}„(*!)&! + (5) 

+ 2! / / {®( qi ~ q2 ">> TT~ u~T A D(t),u)}u(x 1 )u(x 2 )dxidx2, 

where { . , . } is the Poisson bracket. Then taking into account equality (J2J), we observe that the 
distribution functions D n (t), n > 1, are determined by a sequence of the Liouville equations 

^D n (t) = {H n ,D n (t)}, n>\. (6) 

We note that the inverse proposition is also true, namely the generating functional of solutions 
of the Liouville equations (E} satisfies equation (jSJ). By reason of this we are referred to evolution 
equation (jSJ) as the Liouville equation in functional derivatives. 

We consider an equivalent approach to the description of the evolution of states of many- 
particle systems within the framework of a generating functional of marginal distribution func- 
tions [13]. To that end we introduce a generating functional of the marginal distribution 
functions F s (t, xi, . . . , x s ), s > 1, by means of generating functional (jlj) 

(F(t),u) = (D(t),u+l), (7) 

where D(t) = (D(t),I)~ l (l, D,(t, x{), . . . ,D n (t, x±, . . . ,x n ), . . .) is a sequence of normalized 
distribution functions, (D(t),I) is a normalizing factor (grand canonical partition function) 

(D(t), I) = Y]— j ... j D n (t, xi, . . . , £ n )cfei . . . dx n , 

and hence I = (1, . . . , 1, . . .). In consequence of introduced generating functional ([7j) we es- 
tablish the relation between marginal distribution functions and distribution functions. For 
this purpose we transform the generating functional (D(t),u + 1) to form ([T]) to determine the 
sequence of functions generated by such functional. The following equality holds 

(D(t) lU +l) = (e a D(t),u) } (8) 
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where on / e L 1 the operator a (an analog of the annihilation operator [13]) is defined by the 
formula 



(af) n (xx, . . . , x n ) = / f n+x {x 1 ,...,x n+x )dx n+1 . (9) 

Indeed we have 

(D(t),U + 1) = (D(t), I)' 1 J2 — / • • • / D n(t> Xl,..., X n ) Y[(u(x t ) + l)d Xl ...dx n = 

n=0 ■ J J i= i 

oo - « p n n 

= (£>(£), J)" 1 ^— / ... / D n (t,Xi,...,x n )^T ^2 u(x h ) ...u(x ik )dx 1 ...dx n = 

n=0 ' J J k=0 h<-<i k =l 

oo _ oo .. „ „ s 

= (D^),/) -1 ^ — ^ — ■■■ D a+n (t,x 1 ,...,x s+n )Y\u(x i )dx 1 ...dx s+n . 

s=0 ' n=0 ' J •* 1=1 

Therefore, according to (j2j), from ([7]) and (jSJ) we obtain the expression for marginal distribution 
functions within the framework of nonequilibrium grand canonical ensemble 

F s (t,x 1 ,...,x s ) = (D(t),I)- 1 J2^ /••• / Ds+n(t 



3,1, ••• , X s -\- n ) UX s _|_i . . . (IXg-^n. 



Using the Liouville equation in functional derivatives (JSJ), from relation (J7|) we derive the evo- 
lution equation in functional derivatives for the generating functional of marginal distribution 
functions 

j t (F(t),u) = J {K( Pl ), S-WtluMuM + l)d Xl + (10) 

+ 2! j j &* - q ^ 6u( Xl )6u(x 2 ) im ^ n (M(Xi) + l)dXldX2 - 

We are referred to this equation as the BBGKY hierarchy in functional derivatives [6HE]. Tak- 
ing into account equality (J2J), from evolution equation (ITU]) we establish that the marginal 
distribution functions F s (t), s > 1, are determined by the BBGKY hierarchy 



0_ 
Of 



F s (t) = {H s ,F s {t)}+ f {J2${ qi -q s+1 ),F s+1 {t)}dx s+1 , s > 1. (11) 

^ i=i 



As above in case of the the Liouville equations (EJ) the generating functional of solutions of the 
BBGKY hierarchy ( ITTjl satisfies equation ( fTUl ). 

To clarify the structure of stated above evolution equations in functional derivatives we 
generalize them for systems of particles interacting via many-body interaction potentials 

n n n 

H n = J2K( Pl ) + J2 E * ( *W".ftJ> »>1, (12) 

i=l fc=l l=n<...<ifc 
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where $^ is a £;-body interaction potential. In this case the Liouville equation in functional 
derivatives has the form 

±(D(t),u) = J {K(p 1 ),j^-(D(t) } u)}u(x 1 )dx 1 + (13) 

00 Iff s n .". 

+ E~f / ••• / ( $(n) (9i,---,g«.),^^^ -——(D(t),u)}Y[u(x i )dx 1 ...dx nj 

and distribution functions are determined by a sequence of the Liouville equations (E} with the 
Hamiltonian (fT2l) respectively. 

The BBGKY hierarchy in functional derivatives for the generating functional of marginal 
distribution functions takes the form 

j t (F(t),u) = J {K( Pl ), *^F(t),u)}(u(o;i) + 1)^ + (14) 

00 1 r r rn n 

+ g^/--7^'" )te --- 9 " ) - Mx 1 )...Mxj (F(i) '" )} ffi^ +i > fc '--- fc "- 

and the marginal distribution functions are determined by the BBGKY hierarchy respectively 
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F s (t) = {H s ,F s (t)}+j2h E E^/---A $(fe+n)(fe '---' g - ( 15 ) 

fc=l ' i^...^i h =\ n=\ ' J J 

Qs+1, ■ ■ ■ 1 Qs+n), Fs+n{t) }dx s +i . . . dx s + n , S > 1. 



We can see that the structure of the Liouville equation and the BBGKY hierarchy in func- 
tional derivatives is similar in contrast to the structure of these equations for corresponding 
distribution functions. 

2.2 A nonperturbative solution of the BBGKY hierarchy 

On the basis of a solution of the Liouville equation §§§ and relation ([7]) for generating functionals, 
we construct a nonperturbative solution of the Cauchy problem of the BBGKY hierarchy (ITT]) . 
We consider the sequence 

D{t) = (S(-t)D(0), I)-\l, (S 1 (-t)D 1 (0))(x 1 ), ..., {S n (-t)D n (0))( Xl , ..., x n ), . . .), 

where S n (— t), t G R, is the group of evolution operators defined on the space L\ by the formula 

(S n (—t)f n )(xu ...,X n ) = /«(Xi(-t, x u . . . , x n ), . . . , X n (-t, X 1 , . . . , x n )), (16) 

where Xi(—t,xi, . . . ,x n ), 1 < i < n, are solutions of the Cauchy problem of the Hamilton 
equations of n particle system with initial data: Xi,...,x n . The function S n (— t)D n (0) is a 
solution of the Cauchy problem of the Liouville equation §6§ with the initial data D n (0) [13] . 

In the functional (e a S(— t)D(0),u) we expand operators ([TBI) over their cumulants as the 
following cluster expansions 

S 9+n {-t,Y,X\Y)= E IIW-*'^)' n ^ ' ^ 17 ) 
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where Y = (x±, . . . ,x s ), ({Y}) is the set consisting of one element Y = (xi, . . . ,x s ), X = 
(xi, . . . ,x s+n ), and Ep:(m x\y)=U-X- * s ^ ne sum over a ^ possible partitions P of the set 
{{Y}, X \Y) into |P| nonempty mutually disjoint subsets Xi C ({Y}, X \ Y). Owing to 
the equality 

p : ({y},x\r)=u^i **cp zcx\y P:X\y\z=u,i, i ; cp 

and as a result of the symmetry property of the integrand the validity of the following equality 



zcx\y p : x\y\z=u i x i x.cp 



n ?i 



fc=0 ii<...<i fc =l P:(x, +fc+ i,...,x a+ n)=Ui-X'i x ^ cP 



we have 



(^H)D(O), «) = E - E - / . . . / 3W\y|(-t, {F}, I\nEu 



oo ., oo „ „ ,. oo 

1 

X 

,. , , A- 1 

s=0 n=0 J u fe=0 



x J^ Y[^x i \{-t,X i )D s+n+k (0)Y]_u(x i )dx 1 ...dx s+n+k . 

P:(a:«+ T ,+l,...,s s+rl+fc )=Uj.Xj ^CP i=l 

Here the solution of cluster expansions (IT7|) 

2i 1+n (-i,{y}, x\y) = E (- 1 ) |Phl (l p l - 1)« II %^oiH^W) ( 18 ) 

P:({y},x\y)=u i Jf < ^cp 

is the (1 + n)th-oider cumulant of groups of operators f|T6l) . and we introduced the declasteri- 
zation mapping by the following formula: 8({Y}, X \ Y) = X. 
According to the validity of the equality 

/ • • • / E II %\Xi\(-t, Xi)D s+n+k (0)dx s+n+1 . . . dx s+n+k = 

P:(x s + n +l,—,X B +n+k)=\Ji X i X i cP 

= j...JD s+n+k (0)dx s+n+1 ...dx s+n+k , 

and the similar equality for the normalizing factor [13] 

(S(-t)D(0),I) = (D(0),I), 

and taking into account the definition of initial marginal distribution functions, from relation 
(J7J) we derive 

F s (t, Y)=Y2^J-[ a i+-(-*> i Y h X \ W+n(0, X)da; s+1 . . . dx s+n , (19) 



ra=0 
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where Ql 1+n (-t, {Y}, X \ Y) is (1 + n)th-order cumulant (TT8l . 

In fact the following criterion is valid. Expansion ( TT9l) is a solution of the initial-value 
problem of the BBGKY hierarchy (TT5|) if and only if the evolution operators 2li +n (— t), n > 0, 
satisfy the recurrence relations (TT71) . 

We indicate that nonperturbative solution (TT9|) of the BBGKY hierarchy is transformed to 
the form of the perturbation (iteration) series as a result of applying of analogs of the Duhamel 
equation to cumulants (1TB1 of groups of operators [TJ1E3 • 

For the first time a nonperturbative solution of the BBGKY hierarchy was constructed in 
paper [8] (see also [13]) in the form of an expansion over particle clusters which evolution is 
governed by the corresponding-order reduced cumulant of groups of operators f|T6|) 



U 1+n (-t, {Y}, X\Y) = ^(-1) %,, n ,,, & + n-fc(-t, Y, x s+1 , . . . , x s+n . k ). (20) 



fc=0 



k\{n-k)\ 



Indeed, using the reduced cluster expansions of operators (TT6|) over reduced cumulants 



S s+n (-t, Y,X\Y) = ^2 ,,/ J_ ,xj ?7i+fc(-t, {Y}, a; s+x , . . . , a; s+fc ), 

fc=0 '^ '" 

for the generating functional (D(t),u + 1) we have as above 
(D(t),u + 1) = (D(0), I)-\S(-t)D(0),u + 1) = 
= (D(0), I)' 1 J2 j E ^[ / • • • / ^i+n("*) E yD s+n+k (0) H u{x i )dx 1 . . . dx s+n 

s=0 ' n=0 ' ^ ^ fc=0 ' i=l 

and thus, 



F s (t,Y) = f^-f ... I U 1+n (-t, {Y}, Y \ Y)F s+n (0, Y)dx s+1 . . . dx s+n , 

n=0 ' ^ ^ 

where Ui +n (— t) is the (1 + n)th-order reduced cumulant f[2"Uj) . 

3 Evolution equations in functional derivatives for cor- 
relations 

3.1 The Liouville hierarchy for correlation functions 

In addition to an approach of the description of the evolution of states of many-particle sys- 
tems in terms of distribution functions one more an equivalent approach is given by means of 
correlation functions which governed by the Liouville hierarchy. 

We introduce the generating functional (g(t),u) of correlation functions g s (t, xi, . . . , x a ), s > 
1, by means of generating functional (jlj) as follows 

(D(t),u)=e {9{t) ' u) . (21) 
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Then taking into account the equality for k > 1 



Mn) •••Mi*) p :( «,..^) = u^, "p fc(l "'"' fc( Si 

^j = [ x ii i ■ ■ ■ , x i|x • | J 

according to the Liouville equation in functional derivatives (J5]), the evolution of states of many- 
particle systems in terms of generating functional (J2"T1) of correlation functions is governed by 
the following evolution equation in functional derivatives 

i ex 

-(g{t),u) = J {K(p 1 ),—-^(g(t),u)}u(x 1 )dx 1 + (22) 

or in case of the Hamiltonian f[T2"j) this equation has the form 



dt 



d (g(t),u) = J {K( Pl ),j^(g(t),u)}u(x 1 )dx 1 + J2^ J ■■■ J {^ (n) (Qu---,Qn), (23) 

V^ TT <* |Xil (#(*). u) ITT , w 

e n m^.-.m^^ M^ w; 

P:(x 1 ,...,x n ) = U J Xi, **CP V n; V *l^l ; 1=1 

where ^ P .( a . 1 Xn )=U-x- * s ^ ne sum over a ^ P oss ible partitions P of the set (xi, . . . , x n ) into |P| 
nonempty mutually disjoint subsets X,- = (x^, . . . , a?i. x .) C (xi, . . . , x n ). 

Taking into account equality (J2J), we observe that the correlation functions g s (t), s > 1, are 
determined by the Liouville hierarchy [T(3J[TT] 

^(t,F) = {^,^(t,F)}+ (24) 

ipi 

+ E E ■■■ E {^ ¥iliM W--,^pi), n w^)}. 

p : y = U i -X'i, Zicxi, Z| P |CX|p|, x,cp 

|P| > 1 Zi^0 ^|p|^0 

where notations accepted above are used. 

We note that the correlation functions g n (t), n > 1, are interpreted as the functions that 
describe correlations of states D n (t),n > 1, governed by the Liouville equations dHJ). Indeed, 
according to definition fl2]), from relation ( 121]) for generating functionals we obtain the cluster 
expansions of distribution functions over the correlation functions 

D,(t,Y)= Yl liam&Xi), s>l, (25) 

P:Y = {J.Xi XiCP 
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i.e. the correlation functions are the cumulants (semi-invariants) of distribution functions 
9s(t,Y)= J2 (-l) |Phl (|P|-l)! YlDw&Xi), 8>1. 

P:Y = \J i X i XiCP 

The validity of relations ( 125)) is a consequence of the following equality 

eW>'»> = (Exp^t), u), (26) 

where we denote by (Kxp^g(t),u) the generating functional of the sequence of functions 

(Exv,g(t)) a (Y) = J2 U^x^X,), s>l. 

P:Y = \J i X i XiCP 

Indeed, on sequences of functions /, / 6 L 1 we define the ^-product 

(/ * J) m (Y) = J2 f\z\(Z) 7\y\z\(Y \ Z), (27) 

ZcY 

where J2 ZcY is the sum over all subsets Z of the set Y = (xi, . . . , x s ). By means of definition 
(127)) of the ^-product we introduce on sequences / = (0, fi(xi), . . . , f n (xi, . . . , x n ), . . .) the 
mapping Exp^ by the expansions 



oo _. 

Exp,/) |y| (F) = (l + ^-r") |y| (y)= (2f 



n=l 



P:Y=l).Xi X,CP 

where we use the notations accepted above, 8\y\,o is a Kronecker symbol and I = (1, 0, . . . , 0, . . .). 
Then observing the validity of the equality 

(/*/,«) = (f,u)(f,u), (29) 



in view of (128)) we justify equality (Izbj) . 

A solution of the Cauchy problem of the Liouville hierarchy ()24p for the initial data g(0) is 
defined by the following evolution operator [16] 



^(t,y) = 5(t,yb(o))= J2 %i(-*,{*i},...,{*ip|}) U^x^Xi). (30) 

P:V=U,I, XiCP 

Here 2t|p|(— £) is the |P|t/i-order cumulant of groups of operators (fTB")) given by the expansion 
2t| P |(-t,{X 1 },...,{X|p|})= (31) 

e (-lf'^dp'i-i)! n %^i(-*. e ( z *)), 

P':({X 1 },...,{X| P |})=U fc ^ Z fcC P' 
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where Sp' : ({Xi} {x })={] z k * s the sum over an possible partitions P of the set ({X\}, . . . , 
{X|p|}) into |P | nonempty mutually disjoint subsets Zk C ({Xi}, • • • , {X| P |}). 

Really, taking into account formula ( ITEj) for the solution of the Liouville equations and the 
cluster expansions of distribution functions over correlation functions at initial time: -D(O) = 
Exp^g(O), from equality (|26|) we obtain 

(g(t),u) = (Ln*(I + S(-£)Exp^(0)),u), 

where the inverse mapping Ln* to mapping ( 128]) is defined by 

00 ( 1} n ~ * 

(Ln,(i+/)) |y |(y) = (EH - rn \y\ {Y) = 

n=l 

= E (-i) |phl (|p|-i)!ll^iw- 

P:Y=\J t Xi X,CP 

As a result of regrouping terms with respect to the same initial data we derive formula ( 130]) in 
componentwise form. 

3.2 Some relations between generating functionals 

Starting from an alternative approach to the description of states by the generating functional of 
correlation functions (|2TJ) . we define the generating functional of marginal distribution functions 
by means of such generating functional. Moreover we introduce also the generating functional 
of the marginal correlation functions and consider relations of these generating functionals. 

According to relations ([7|) and (T2T1) . the generating functional of marginal distribution func- 
tions is determined by means of generating functional of correlation functions as follows 

(F(t),u) = e {9{t) ' u+1) - ( - 9{t)J) . (32) 

Using this relation, the BBGKY hierarchy in functional derivatives (TT4l) can be derived on the 
basis of the Liouville hierarchy in functional derivatives (1221) . 

From (1321) in view of definition (J2J) we obtain the expression for marginal distribution func- 
tions in terms of correlation functions 

F s (t,Y) = y^~ / • • • / 9i+n(t, {Y}, x s+1 , ..., x s+n )dx s+1 . . . dx s+n , s > 1, (33) 

where we denote by ({^}) the set consisting of one element Y — (1, . . . , s), and the correlation 
functions gi+ n (t, {Y},x a+ i, . . . ,x s+n ), n > 0, are solutions of the Liouville hierarchy ffM§ in 
case of particle cluster of s particles and particles [18]. To prove relation (|33|) we transform the 
functional from the right-hand side of relation (|32l) to the canonical form using equalities ( 126|) 
and (JHD 

e G,(t),«+i)-to(t) 1 i) = (Exp^(t),/)- 1 (Exp^(t), M + l) = 
= (Exp,g(t),I)- l (e a Ex V ,g(t),u). 
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The last functional is the generating functional of the sequence of functions 

(Exp^(£),/)-VExp^(£). 

To write down it in the component- wise form on measurable functions we introduce the following 
mappings 

(^Yf)n — f\Y\+n{Y, X s +i, . . . , X s+n ), 

(Q{Y}f)n = /l4«({Y}, Xs+i, . . . , X s + n ), Tl > 0. 

Then we have 

(Exp^(t),/)- 1 (e a Exp^(t)) s (F) = (Exp^(t),/)- 1 (O y Exp^(t),J). 
Owing the validity of the following equalities [IB] 

OyExp^(t) = D {y} Exp„0(£), 
and 

n {y} Exp^(t) = Exp^(t) * d {Y }g(t), 
and according to equality (1251) . we finally derive (135)) 

(Exp^(t),/)- 1 (O y Exp^(t),J) = (d {Y} g(t),I). 

On the basis of solution (l30l) of the Liouville hierarchy (1221) and relation (1321) for generating 
functionals, we also derive nonperturbative solution (1T91) of initial- value problem of the BBGKY 
hierarchy (1TT1) . 

By analogy with relation (j7j) we introduce the generating functional of the marginal corre- 
lation functions by means of generating functional of correlation functions as follows 

(G(t),u) = (g(t),u + l) - (g{t),I), (34) 

or in view of definition ([2]) in componentwise form it has the form 

G s (£, Xi, . . . , x s ) = y — I ... I g s+n {t, X\, . . . , x s+n )dx s+ i . . . dx s + n , s > 1, (35) 
n=0 ' ^ ^ 

where the correlation functions g s+n (t) are solutions of the Liouville hierarchy (124"]) . We em- 
phasize that every term of marginal correlation function expansion (1351) is determined by the 
(s + n)t/i-particle correlation function as contrasted to marginal distribution function expansion 
(13"3"|) which is defined by the (1 + n)£/i-particle correlation function. 

From relations (134"]) and (132]) we obtain the relation of the generating functionals of marginal 
distribution functions and marginal correlation functions 

(F(t), M ) = e (G(t) ' u) . (36) 

We note that in view of definition (J2J) in componentwise form relation ( 136]) is the cluster ex- 
pansions of the marginal distribution functions over the marginal correlation functions (it is an 
analog of cluster expansions ( 125]) or ( 12T]) ) 

F s (t,Y)= J2 II G l*l(*' X *)' s ^ lj 

P.Y = \J i X i X lC P 

where Y1p-y=\j. x t * s ^ ne sum over a ^ possible partitions P of the set Y = (xi, . . . ,x s ) into |P| 
nonempty mutually disjoint subsets Xi C Y. 
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3.3 The nonlinear BBGKY hierarchy in functional derivatives 

From relation (|34"|) we derive the nonlinear BBGKY hierarchy in functional derivatives for the 
generating functional of marginal correlation functions on the basis of the Liouville hierarchy 
in functional derivatives (1221) 



j t (G(t),u) = J {K(px), J A-}(G(t),u)(u(x 1 ) + l)d Xl + (37) 

or in case of the Hamiltonian f|T2|) this equation has the form 



J n C 

-(G(t),u) = J {K( Pl ),j—r(G(t),u)}(u(x 1 ) + l)dx 1 + (3* 

n=l J ^ P:f.x, ^ = II.X,. XvCP°^ X * 1 ^-- ^ a; W 






n 
X 

4 = 1 



PJ(«(a?i) + l)dxi. . .dx n , 



where we use notations accepted above in equation (1231) . 

Taking into account equality (J2J), we observe that the marginal correlation functions are 
determined by the nonlinear BBGKY hierarchy 

^G,(t,Y) = {H„G.(t,Y)} + (39) 

+ E E ••• E {^ |Zrl) (^i,...^iPi), n<W'*)} + 

p : y = U l ^ l , ZicXi, Z|p|CX| P |, x.cp 

|P| > 1 Zi^0 Zipi^O 



4-n+l-fcJ, 



+EE E [...[{® {n+1) (<in,...,<i Jk ,<is + i,..., qt 

n=l fc=l l=ji<...<j fc ^ ^ 

P:(xi,...,x s+n+1 _ fc ) = U J X J ,|P| <n + l, X.CP 

x i ?r\( % ,...,i js ) 

We can see that the structure of the Liouville hierarchy ([22]) and the nonlinear BBGKY 
hierarchy (|37|) in functional derivatives is similar in contrast to the structure of these hierarchies 
for corresponding correlation functions. 

On the basis of solution (1301) of the Liouville hierarchy (1241) and relation (1341) for generating 
functionals, as above in case of relation (J7J) we derive a nonperturbative solution of initial- value 
problem of the nonlinear BBGKY hierarchy (139]) . 
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Taking into account the equality 

(g(t), u + 1) = (S(t\g(0)),u + 1) = (e a S(t\g(0)),u), (40) 

in the functional (e a S(t\g(0)),u) we expand operators ( 1301) over the corresponding cumulants. 
As a result from relation ( 1341) we derive a solution expansion of the nonlinear BBGKY hierarchy 

m 

G s (t,Y) = y^~ / • • • / ^i+n(*5 {^}, £ s +i, • • • , a^+n | G{0))dx s+1 . . . dx s+n , S > 1, (41) 



n=0 



where the (n + l)£/i-order reduced cumulant Ui +n (t) of groups ( 130]) has been introduced 
U 1+n (t; {Y}, x s+1: ..., x s+n | G(0)) = X)(-l)* , " , x (42) 



fc=0 



fc!(n-A;)! 



P:(xi,...,x s+n _ fc )=|J^i fc l=0 fe| P | 

X^IXil+fc-fcilU) ^-li x s+n-k+li ■ ■ ■ i x s+n-ki) ■ ■ ■ G |X P i |+fc|p|_i (0, X|p|, ^s+n-tpii+b • • ■ ; ^s+n) 



— ' fci!(fc - fei)! ' ^ fc|p|_i!(fc|p|_ 2 - fc|p|-i) 

1=U fc,p, 1 =U 



and 2l|p|(— t) is the |P|£/i-order cumulant defined by the formula (l3Tj) . 
We give simplest examples of reduced nonlinear cumulants 



U 1 {t;{Y}\G(0)) = g(t;Y\G(0)) = 

= Yl *\*\(-t,{Xih'-;{*\p\}) II <W°.*0. 

U 2 (t;{Y},x s+1 \G(0))= Yl 2t| P |H,{^i}, ■■-,{*|P|}) II <W°>*<)- 

P:(y,a!, + i)=Ui-X'i XiCP 

|P| 

£ %,(-*, {*i},---,{*|p|})£ J] G w (0 ) X i )G|^| +1 (0,X i) x H . 1 ). 
P:y=U l x< j=i x, c p, 

Xi ^t Xj 

In case of initial data satisfying a chaos property, i.e. G^(0) = (0, Gi(0, 1), 0, . . .), for the 
^1 + n)th-OTder reduced cumulant (J4"2l) we have 



*7 1+n (t; {Y}, X\Y\ G (1) (0)) = 2l s+n (-t, X) 11 Gi(0 






i=l 



where 2l s+n (— t) is (s + n)t/t-order cumulant f|T8l) . and solution expansion 04 ip takes the more 
transparent structure. From the structure of this series it is clear that in case of absence of 
correlations at initial time the correlations generated by the dynamics of many-particle systems 
are completely governed by cumulants ( TT8l) of groups of operators ( ITBl . 
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4 The dual BBGKY hierarchy in functional derivatives 

4.1 A generating functional of marginal observables 

The evolution of many-particle systems is described within the framework of the evolution 
of states or in terms of an equivalent approach as the evolution of observables [13]. The 
evolution of observables of a system of a non-fixed number of particles is described by the 
sequences A(t) = (A , Ax(t, Xi), . . . , A n (t, x%, . . . , x n ), . . .) of observables. The bounded function 
A n (t, Xi, . . . ,x n ) is defined on the n-particle phase space and symmetric with respect to the 
permutations of arguments xi, . . . ,x s . According to definition ([1]) we introduce the generating 
functional (A(t), u) of observables by the expansion 

{A(t),u) = ^2~ ■■■ A n (t,x 1 ,...,x n )Y\_u(x i )dx 1 ...dx n . (43) 

For the integrable functions u(xi), i > 1, functional ( I43p exists. 

Within the framework of generating functional ( 143]) the evolution of observables of many- 
particle systems with the Hamiltonian (]3]) is governed by the following evolution equation in 
functional derivatives 



1 p C 

j t {A{t),u) = J {——(A(t), u), K( Pl )}u(x l )dx 1 + 

+ 2! / / { fT \ S ( d A(t),u)^(q 1 -q2)}u(x 1 )u(x 2 )dx 1 dx2, 



(44) 



where { . , . } is the Poisson bracket. Taking into account equality (j2]), we establish that observ- 
ables A n (t), n > 1, are determined by the sequence of the Liouville equations for observables 

^A n (t) = {A n (t),H n }, n>l, (45) 

or vice versa the generating functional of solutions of the Liouville equations (I45p satisfies 
equation (T44]) . 

We introduce a generating functional of marginal observables B(t) = (B Q , Bi(t, Xi), . . . , 
B s (t, Xi, . . . , x a ), . . .) by means of generating functional (143]) as follows 

(B(t),u) = (A(t),ue-^^ d n, (46) 

where we use the notation 

ue -fu {x)dx = ( i ; m(xi); . . . ? Y[ u (xi), ...)J2 -jr ■■ U< x ^ dXi ■ ■ ■ dXk - 

i=l fc=0 ' J J j=l 

A generating functional of marginal observables gives an equivalent approach to the description 
of the evolution of finitely many particles in addition to generating functional of observables 
and it is adopted for the description of infinite-particle systems. 
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In consequence of introduced relation ( 146J) we establish the relation between marginal ob- 
servables and observables. We transform generating functional (A(t),ue~-> u ( x > dx } to form (CQ), 
and thus, we determine functions generated by this functional. The following equality holds 
(the dual analog of equality ([8])) 

(A{t),ue-M* )dx ) = {e- a+ A{t),u), (47) 

where on b G L°° the operator o + is defined by the formula (an analog of the creation operator 
|13j . which is an adjoint operator to operator (Q in the sense of mean value functional [13] ) 

s 

(a + b) s (xi,...,x s ) = 22 b s -i((xi,...,x s ) \ (xj)). 

3=1 

Indeed we have 

oo _. oo / ^j, „ „ k+n 

(A(t),Ue-f u(x)dX )=J2^jJ2kT / " A nit^U---^n)\[u(x i )dx 1 ...dx k+n = 
n=0 ' k=0 ■ J J i=l 

oo _. „ „ n ■ n 

Yl ~\ / • • • / Y^(~ l } k k\(n- k\\ An - k ^ Xl ' " ' ' Xn ^ II u ^) dx i ■ ■ ■ dx ^ 



n=0 " " fc=0 v ' i=\ 



and therefore in view of symmetry property of observables, according to definition (|2j), from 
relation ( H6l) we obtain the expression for marginal (s-particle) observables within the framework 
of nonequilibrium grand canonical ensemble [20] 



s ( _-\\n s 

B a (t,x l ,...,x a ) = '^2 — $^ A B-n{t,(x 1 ,...,x„)\(x jl ,...,x jn )). (48) 

Using relation (146 p . from the Liouville equation in functional derivatives for observables 
we derive the evolution equation for the generating functional of marginal observables 



B(t),u) = J {jJ—-(B(t),u) J K(p 1 )}u(x l )dx 1 + (49) 



C2 2 



+ UI l M^kw {m - m) + £ fafa {m u)Mqi - 92)} n -w^, 



8=1 



We are referred to this evolution equation as the dual BBGKY hierarchy in functional deriva- 
tives. Taking into account equality ([2]), from evolution equation (T4T)]) we establish that the 
marginal observables B s (t), s > 1, are determined by the dual BBGKY hierarchy 



Ft s 

B s (t,Y) = {B s (t,Y),H s }+ J2 {B-i&YMxjmqi-qj)}, s > 1, (50) 



at 



where Y = (xi, . . . , x s ) and F \ (xj) = (xi . . . , Xj_i, Xj+i, . . . , x s ). As above in case of the 
Liouville equations (]45p the generating functional of solutions of the dual BBGKY hierarchy 
(150]) satisfies equation 
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To clarify the structure of stated above evolution equations in functional derivatives we 
generalize them for a system of particles interacting via many-body interaction potentials. In 
this case the Liouville equation in functional derivatives has the form 

f t {A(t),u) = J {——{A{t),u),K{p 1 )}u{x 1 )dx 1 + (51) 

00 1 r r rn n 

+ g^y ■■■] ^ ^..i^ (i(f ^^ (n)fe ----' 9 ^n^^-'-^ 

and observables are determined by a sequence of the Liouville equations ( l4"5j) with the Hamilto- 
nian f|T2l) respectively. The dual BBGKY hierarchy in functional derivatives for the generating 
functional of marginal observables takes the form 

-(B(t),u)= {——(B(t),u),K( Pl )}u(x 1 )dx 1 + J2^ x ( 52 ) 

/„ n n <-fc n 

... {J2 E J^U~\ T-7-T(B(t)^),^ n) (q 1 ,...,qn)}l[u(x i )dx 1 ...dx n , 
J k=l h<...<i k = l V 'I/ --- V l fc/ i=l 

and the marginal observables are governed by the dual BBGKY hierarchy [2T] 

^B s (t,Y) = {B s (t 1 Y) 1 H s }+ (53) 

s , s _. s 

+ E^ E 7fcT^ S {^-n(t,r\(x J1 ,...,X,J),$«(x J1 ,...,X J j}. 

As is obvious the evolution equations in functional derivatives for dual generating functionals 
P|). p3j) have the form (fT3l . (l5"Tl) and correspondingly for dual generating functionals fl7|).( |46l) 
the hierarchies in functional derivatives have the form flT4|) . fl52|) . 

4.2 A nonperturbative solution of the dual BBGKY hierarchy 

On the basis of a solution of the Liouville equation ( H5j) and relation (T46T) for generating func- 
tionals, we construct a nonperturbative solution of the Cauchy problem of the dual BBGKY 
hierarchy (|53|) . With this aim we consider the sequence 

A(t) = (Ao, S l (t)A 1 (0, xi), . . . , S n (t)A n (0, Xl , . . . , x n ), . . .), 

where the evolution operator S n (t) is defined on the space L£° by formula (TTOl) . and it is 
determined a solution of the Cauchy problem of the Liouville equation for observables of 
the n-particle system with the initial data A n (0). Using equality (JSj), in the functional 
(e~° S(t)A(0), u) we expand operators S n (t), n > 1, over their cumulants as the following 
cluster expansions 

(-l) n S s . n {t,Y\X) = (54) 

n 

= E( _1 ) n ~ fc E %i+k(t,{Y\(x il ,...,x ih )},x il ,...,x ik ), 

fc=o n<...<j fe e(ii,...j„) 
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where the abridged notations are introduced: X = (xj 17 . . . , Xj n ) C Y, (x^, . . . , x ik ) C X. Then 
the following equalities are true 

(e- a+ S(t)A(0),u) = 

= J2 1 / • • • / E ^T ^ 5 s _„(0A_ n (0)(y \ X) f[ u(x i )dx 1 ...dx s = 

s=0 " J J n=0 " j 1 ^...^ n =l 

x ^ ^ A,_„_ fc (0,F\X\ (a^,...,^)) JJ«(a;i)ctei...cte s . 

fc=o ' ^^...^^ = 1, i=l 

111-")** 7^ ill- ■■:J , n 

Owing to the definition of initial marginal observables ( I48p . from relation (I46p and equality 
( 14 7p we finally derive 

S -. s 

B s (t,Y) = J2- } E 2li+n(t,{y\X},X)i? s _ n (0,r\X), (55) 

where X = (xj 1 , . . . , x Jn ) C K and the (1 + n)th-order cumulant of groups of operators ( JTTjj) is 
defined as follows [15], EH 



2l 1+n (t,{F\X},X)= £ (_i)|P|-i(| P |_i)! II ^)iM(X)), (56) 

p : ({y\x},x)=u^ X lC P 

and we use notations accepted above, for example, 

Mt,{y}) = s s (t,Y), 

2l 2 (t,{F\ (x,)},x,) = S s (t,Y) - 5 s _!(i,y \ (xj^S^Xj). 

To prove that cumulants (156]) are solutions of recurrence relations ( l54"j) we observe that 
recurrence relations ( )54|) hold on functions b s - n (Y\X), and thus, the following identity is true 

(-l) n S s . n (t,Y\X) Y, (- 1 )' P| I P I ! II S lXil (t,X t )b s ^(Y\X) = 

P:X=\JiXi X,CP 

n 

= Y^(-l) n - k Y, %i +k (t,{Y\(x h ,...,x ik )},x il ,...,x ik )b s „ n (Y\X), 

k=0 u<...<i fe G(ji,...,j„) 

where Ylpx=\\x * s the sum over an possible partitions P of the set X into |P| nonempty 
mutually disjoint subsets Xi C X. These recurrence relations are rearranged form of the 
cluster expansions of the evolution operators S s (t), s > 1, 

S s (t,Y\X,X)= J2 n%'l(*' X *)' l X l- ' ( 57 ) 

P:({Y\X},X)=U i X i X Z CP 
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and hence we obtain expression ( 1561) as a solution of recurrence relations ( 1571) . 

In fact as above the following criterion is valid. Expansion ( 15*5*1) is a solution of the initial- value 
problem of the dual BBGKY hierarchy ( 153]) if and only if the evolution operators 2li +n (£), n > 0, 
satisfy the recurrence relations ( 157|) . 

We remark that an analog of expansion (TT9|) - (|20|) in case of a solution of the dual BBGKY 
hierarchy (15T)]) is given by expansion (1551) with the following reduced cumulants of groups of 
operators (""J"""]) 

n | 

U l+n (t, {xx, . . . , x s _ n }, x s _ n+ i, ...,x s ) = ^2(- 1 ) k kl(n ^ k)] s s-k(t). (58) 

fc=0 '^ >' 

In fact it is a consequence of the equivalence of corresponding generating functionals. Indeed, 
solving recurrence relations (15TJ) with respect to the first-order cumulants for the separation 
terms which are independent from the variables Y \ X 

p 

a 1+n (t, {y \ x}, x) = J2 Mt, {y\xuz}) J2 (" 1 ) |P| l p l ! Tl^ (*> { x *})' 

where J2zcx * s ^ ne sum over an P oss ible subsets Zclof the set X, and, taking into account 
the identity 

IP 
Y, ("1) |P| |P|! H*&> {Xi})b^n{Y \X)= J2 (" 1 ) |P| \ P \ ] -b s - n (Y \ X), 

P:X\Z=(J t Xi i=l P:X\Z=\J.Xi 



and the equality 



£ (-1)' P '|P|! = (-l)'^l 



we derive expansion (153]) over reduced cumulants (1581) of the dual BBGKY hierarchy (1501) . 

We note that for certain classes of observables solution expansion (1551) has more simple struc- 
ture. As we can see from (|48j) one component sequences of marginal observables correspond to 
observables of certain structure, namely the marginal observable B^ = (0, ai(x\), 0, . . .) corre- 
sponds to the additive-type observable A^ = (0, a 1 (s 1 ), . . . , Y^=i a ii. x i)i ■ ■ ■)■> an d in the general 
case the /c-ary-type marginal observable B^ = (0, . . . , 0, a k (xi, . . . , Xk), 0, . . .) corresponds to 
the fc-ary-type observable A^> = (0, . . . , 0, a fc (xi, . . . , x k ), . . . , J27 x <...<i k =i a s(x h , . . . ,x ik) , . . .)• 
Then in case of additive-type observables solution expansion (155]) attains the form 

s 

Bi 1 \t,Y)=Ql s (t,Y)J2a 1 (x l ), 

i=l 

where 2l s (£, Y) is the sth-order cumulant (156]) of groups of operators ("It"") . 
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5 Conclusion 

By means of established relations for generating functionals (IT]). (l2~T]) .( l3"4"]) and fj46|) the hier- 
archies of evolution equations in functional derivatives for generating functionals of marginal 
states and observables of classical many-particle systems, namely the BBGKY hierarchy for 
marginal distribution functions ( ITUj) . the Liouville hierarchy ( 122]) for correlation functions, the 
nonlinear BBGKY hierarchy for marginal correlation functions ( )37|) and the dual BBGKY hi- 
erarchy for marginal observables (1521 . were derived from the Liouville equations in functional 
derivatives for generating functionals of distribution functions (]5]) and observables (joTj) respec- 
tively. To clarify the structure of these hierarchies in functional derivatives we generalized them 
for particles interacting via many body interaction potentials, in particular we have equations 
(ir3"]) .( ]T4"]) for states, equations (I2"3"]) .( 13"5]) for correlations of states and equations (15"1~]) . (15"2"]) for 
observables. 

On the basis of cluster expansions ffTTl) and (154]) of groups of operators which describe 
dynamics of finitely many particles the nonperturbative solutions of the hierarchies of evolution 
equations for states and observables were constructed. They are represented in the form of 
expansions ( TTT?|) and ([57)]) over particle clusters which evolution is governed by the corresponding 
order cumulant (TT8|) of groups of operators (TT6|) of the Liouville equations (Q in case of the 
BBGKY hierarchy ( 1T5J1 and cumulant (|56|) of groups of operators of the Liouville equations for 
observables (]4"5]) in case of the dual BBGKY hierarchy ( 153]) respectively. 

Moreover, along with the definition within the framework of the non-equilibrium grand 
canonical ensemble the marginal distribution functions can be defined in terms of dynamics of 
correlations ( 132]) that allows, in particular to give the rigorous meaning of the states for more 
general classes of functions than the integrable functions. Using relation ( ]34"]) of the generat- 
ing functional of marginal correlation functions and the generating functional of correlation 
functions governed by the Liouville hierarchy f )24|) . a nonperturbative solution of the nonlin- 
ear BBGKY hierarchy for marginal correlation functions (13"9"]) was constructed in the form of 
expansions (J4T]) over particle clusters which evolution is governed by the corresponding order 
cumulant ( 14~2]) of groups of operators (130]) of the Liouville hierarchy ( 124]) . 

Thus, the concept of cumulants of groups of operators forms the basis of established nonper- 
turbative solution expansions of hierarchies of evolution equations of many-particle systems. 
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